Superfluid— Insulator Transition in Commensurate One-Dimensional Bosonic System 

with Off-Diagonal Disorder 
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We study the nature of the superfluid-insulator quantum phase transition in a one-dimensional 
system of lattice bosons with off-diagonal disorder in the limit of large integer filling factor. Monte 
Carlo simulations of two strongly disordered models show that the universality class of the transition 
in question is the same as that of the superfluid-Mott-insulator transition in a pure system. This 
result can be explained by disorder self-averaging in the superfluid phase and applicability of the 
standard quantum hydrodynamic action. We also formulate the necessary conditions which should 
be satisfied by the stong-randomness universality class, if one exists. 
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An interplay between commensurability, interactions, 
and disorder in the superfluid-insulator quantum phase 
transition is a challenging theoretical problem with ap- 
plications to such diverse physical systems as ^He in 
porous media and aerogels, superfluid fllms on various 
substrates, Josephson junction arrays, granular super- 
conductors, disordered magnets, etc. (see, e.g., P, 0, 
0,0, [StlSlJ, and references therein). Recently, research 
on disordered bosons is strongly stimulated by the pos- 
sibility of controlled experiments with ultracold atomic 
systems |Hi|- 

Commensurability is relevant for the system of lat- 
tice bosons with the off-diagonal disorder (random hop- 
ping amplitude or on-site repulsion) and exact particle- 
hole symmetry (taking place in the limit of large fill- 
ing factor) 1^]. It leads to a new type of insulator, in- 
compressible and gapless Mott glass (MG) Ji)J in place 
of incompressible gapped Mott insulator (MI). In two- 
dimensional (2D) systems, changing the nature of the in- 
sulating phase results in the new universality class of the 
superfluid-insulator (SF-I) transition. The SF-MG tran- 
sition is characterized by the dynamical critical exponent 
1 < z <2 different from z = 1 for the SF-MI point in a 
perfect system 0]. A similar situation is expected in 3D. 

Surprisingly, very little is known for fact about the 
ID case apart from perturbative renormalization group 
(RG) arguments. On the superfluid side of the transi- 
tion, one can use the instanton language in terms 
of which weak off-diagonal disorder does not seem to be 
relevant, leading only to the inhomogeneity of the mi- 
croscopic stiffness in the equivalent (H-l)D classical XY 
model. Recently, Altman et al. argued on the basis of 
the spatial RG analysis that while small off-diagonal dis- 
order is indeed irrelevant for the SF-MG criticality, there 
should exist also the strong-randomness fixed point. We 
are not aware of any large-scale numerical simulations at- 
tempting to study the SF-MG criticality in the strongly 
disordered system. 



The main observable of interest in ID systems is the 
Luttinger-liquid parameter g = ttVAk, where A is the 
superfluid stiffness and k is the compressibility. There is 
a "smoking gun" signature in the behavior of g at the 
critical point which allows to discriminate between dif- 
ferent scenarios of the phase transition. If the SF-MG 
transition is in the same universality class — at least on 
the superfluid side — as the SF-MI point in a pure sys- 
tem then g should jump discontinuously from gc = 2 
to zero and obey the Kosterlitz-Thouless RG equations 
in finite-size systems (the corresponding jumps in A and 
K are not universal). The prediction of Ref. 6] for the 
strong-randomness fixed point is that only k obeys the 
Kosterlitz-Thouless (KT) scenario Our approach 

then is to simulate the (1+1)D classical analog of the 
bosonic Hubbard Hamiltonian using Worm Algorithm 
and to see whether the KT scenario with the uni- 
versal jump = 2 is taking place. 

We simulate two microscopically different models of 
strong off-diagonal disorder. In both cases we observe 
agreement with the ideal KT scenario. This fact and 
the observation that A and k are self-averaging quanti- 
ties near the transition point allow us to put forward an 
asymptotic quantum hydrodynamic approach. We argue 
that instantons of unit charge inevitably become relevant 
when g is small enough and thus any consistent strong- 
randomness scenario should necessarily predict discon- 
tinuous behavior of all three quantities A, k, and g. 

At large integer filling factor, the bosonic Hubbard 
model is equivalent to the quantum rotor Hamiltonian 
(see, e.g., [i3|): 



H = -J2t, cos(^,+i - ^^.) + E T (7 9^ 



(1) 



The first term describes particle hopping between the 
nearest neighbour sites and the second term describes 
the on-site repulsion. 




FIG. 1: Luttinger-liquid parameter g for the bimodal distri- 
bution of tx- Statistical errors are smaller than the symbol 
size. The solid line is the KT extrapolation to the infinite 
system size. Dotted lines are to guide an eye. 



While it is possible to perform a direct simulation of 
the quantum model we prefer to work with its macro- 
scopic classical equivalent — a version of the bond-current 
model 01 described by the anisotropic action: 



2 ^ 



(2) 



Here integer vector n — (x, r) labels sites of the two- 
dimensional square space-time lattice of linear size L, x 
and f are unit vectors pointing in the space and discrete 
imaginary time directions, respectively, and Jn,Q {a = 
x,f) are integer "currents" living on lattice bonds. The 
allowed configurations of bond currents are subject to the 
zero-divergence constraint X]a('^n,c( + </n,-a) = 0, where, 
by definition, the direction —a is understood as opposite 

to a and Jn.-a = —Jn-a,a- 

As its limiting case, the bond-current model |(2Jl repro- 
duces the (1+1)D worldline action of the Hamiltonian 
the correspondence being given by the relations 
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U = U-At , 



(3) 



where At is the imaginary time step and the limit of 
At — *■ is assumed. The trick, however, is to consider 
essentially finite At; this may not change the universality 
class of the continuous phase transition but substantially 
improves the algorithm performance 0, 0, 0| ■ 

For each system size we considered 2 • 10^ -^2-10'^ 
realizations of disorder. The analysis of sample-to-sample 
fluctuations shows that there is self-averaging of A and k 
near the SF-MG transition for both models of disorder. 
Qualitatively, the behavior of A and k is similar to that 
of g and indicates the presence of a jump at the critical 



FIG. 2: RG parameter ^ for different data sets as a function of 
control parameter U for the bimodal distribution of tx- The 
solid line is a linear fit. 



point. In what follows, we discuss only the Luttinger- 
liquid parameter g — tt^/Kk because its jump is expected 
to be universal and can be deduced from the quantitative 
RG analysis of finite-size effects. 

The first set of simulations was performed for fixed 
Ux = U (i.e. no disorder in the on-site interaction) and 
the bimodal distribution for the hopping parameter t^'. 
with equal probabilities we choose = 21} or tx — 2C//3. 
Simulation results and comparison with the RG flow for 
g as, a. function of system size L are presented in Figs.^ 
121 When reduced to its integral form, the KT renormal- 
ization flow is given by 



ff(il)/2 



dt 



(L.)/2 t^{\nt-^) + t 
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(4) 



where ^ is size-independent microscopic parameter char- 
acterizing the vortex fugacity. The procedure of ana- 
lyzing the data is as follows. For different pairs of sys- 
tem sizes, Li and L2, and fixed value of U close to the 
critical point, one solves Eq. for ^, and verifies that 
for large system sizes ^ is L-independent. Moreover, the 
S,{U) dependence must be analytic and thus well approx- 
imated by a straight line in the vicinity of the critical 
point U = Uc- In contrast, the gL{U) curves noticeably 
deviate from the straight line even in a rather small in- 
terval around Uc, as they develop a square-root cusp and 
the universal jump at C/ — Uc in the limit of i 00. 
These features are clearly seen in Figs. ^ |21 

It seems that the bimodal distribution with the ra- 
tio of hopping amplitudes as large as three may be re- 
garded as a strong disorder case since our system sizes 
accomodate many local disorder fluctuations with 3 5 
consequtive bonds with small/large t. Still, one may not 
exclude the possibility that weak links with power law, 
rather then exponential, distribution function of small t 




FIG. 3: Luttinger-liquid parameter g for the power-law dis- 
tributed tj. Statistical errors are smaller than the symbol 
size. The solid line is the KT extrapolation to the infinite 
system size. Dotted lines are to guide an eye. 



is necessary to realize the strong-randomness fixed point. 
We thus have studied a system with continuous distribu- 
tion of links between the nearest neighbour sites. The 
quantum-rotor prototype for our model was the Hamil- 
tonian Q with the j-independent Uj = U, and hop- 
ping amplitudes tj independently distributed on the in- 
terval (0, ^rnax) with the power-law probability density 
p{tj) (X ij, To emphasize a possible special role 

of having more frequent weak links, we kept the ratio 
tmani/U constant and used 7 as a control parameter. The 
closer is 7 to zero, the stronger is the effect of weak 
links. Accordingly, for the bond-current model ((SJ we set 
Ux = 0.5 and tx — —2 ln(C/2), where C is a random num- 
ber distributed with the probability density p(Q oc on 
the (0,1.5) interval. 

In Fig. 121 we present our data for g as a function of 7, 
which is qualitatively similar to that in Fig.Q] Note that 
strong finite-size effects set in only at < 2. Also, we 
observe that self-averaging of g requires very large system 
sizes L > 100 in the vicinity of the transition point. This 
explains why fitting the data with Eq. Q), see Fig. 01 
demonstrates stronger finite-size corrections than in the 
previous case. The overall behavior is consistent with the 
pure-system KT scaling. 

Self-averaging of mesoscopic A and k was observed in 
both simulations close to transition point for g > 2. On 
the superfluid side of the transition, this outcome is ex- 
pected and can be proved under two quite general as- 
sumptions. If one takes two ID systems of size L (label 
them as systems 1 and 2) and combines them together 
to form a new system of size 2L, then k{2L) and A{2L) 



FIG. 4: RG parameter ^ for different pairs of data sets as a 
function of 7 for the power-law distributed tj. The solid line 
is a linear fit. 

are given by 

At(2L) = [ki(L) + K2(L)]/2, (5) 
A-'{2L) = [A^\L) + A^\L)]/2. (6) 

Equations ((SJ and 10 imply that properties of the new 
system are independent of the junction properties. The 
latter assumption can hardly be questioned if initial k's 
and A's are finite and L is large enough. The self- 
averaging of K is then guaranteed by the Central Limit 
Theorem (CLT). The case of A is less obvious, since the 
distribution function f{L,A) may be such that (A~^) = 
00, in which case the CTL does not apply. However, if 
the infinite-size limit for the distribution function is well 
defined, f{L 00, A) = /(A), then it must satisfy the 
following relation 

/(A) = / <5 (a - ^^^) /(Ai)/(A2) dAidA2 . (7) 

Now it is easy to prove that the only solution of Eq. 
is the self-averaging one, /(A) = S{A — Aq). Indeed, 
Eq. O leads to the inequahty (A^) < (A)^, immediately 
implying that the distribution is dispersionless. 

When the self-averaging of A and k takes place, the 
long-range groundstate properties of the superfluid phase 
can be described by the Popov's hydrodynamic action 
[itI i with the density-field fiuctuations integrated out: 

Sm= |dxdr(^^{<^>',r + ^iKf+^no(.x)K]■ (8) 

Here ^{x,t) is the phase field, and uq is the equilibrium 
value of the local number density. The last term in S 
is sensitive only to instantons — vortex-type topological 
defects in the non-single- valued phase field $(a;,r). In- 
tegrating the last term over r and introducing a rescaled 
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time variable y = ct {c = a/A/k is the sound velocity), 
one gets the effective action in the form (see, e.g., 

Sm = ^J {^<^fdxdy + iY,p, e{x,) , (9) 

where v enumerates the mstantons, = ±1, ±2, . . . and 
are the charge and the coordinate of the instanton v, 
respectively, and 

9{x) = -2tt I na{x')dx' . (10) 
Jo 

The first term in Eq. ^ corresponds to the 2D XY 
model, and the second term describes the position- 
dependent instanton phase. For a system in a periodic 
external potential, one can ignore the phase term if it 
changes by a multiple of 27r when the instanton is shifted 
by one lattice period, a, i.e. when for any integer m 

/ no{x)dx = integer. (11) 
Jo 

This relation is satisfied for a pure system at integer 
filling factor, but typically does not hold true in the 
presence of disorder because no{x) becomes a random 
function of x. An exceptional situation occurs in a sys- 
tem with off-diagonal disorder and at large integer fill- 
ing. The particlc-holc symmetry is then preserved locally, 
and the occupation number for each site x remains inte- 
ger. Hence, the second term in @ is irrelevant and the 
macroscopic behavior of the system is indistinguishable 
from that of a pure system, i.e. the SF-MG transition is 
in the 2D XY universality class. 

Under these circumstances, one can make the follow- 
ing statement about alternative (to the pure-system one) 
scenario of the SF-MG criticality. It is possible only if at 
the critical point it leads to a jump of g with the ampli- 
tude greater or equal than 2. Moreover, at gc the system 
should be non-self-averaging, to exclude the applicability 
of the hydrodynamic approach. Any theory that predicts 
a SF phase with finite A and k and g < 2 is inconsis- 
tent, since, in view of the established self-averaging, the 
hydrodynamic+instanton approach is applicable and su- 
perfluidity is destroyed by the KT scenario at g = 2. 

It is worth noting that while the hydrodynamic action 
© is applicable to the SF phase and the critical point, it 
can not be used on the insulating side of the transition. 
Formally, for g < 2 Eq. © predicts the MI phase with 
the gap in the energy spectrum while the MG phase is 
gapless (the absence of the gap is due to arbitrary large 
and thus exponentially rare superfluid regions |^). 

To conclude, we have presented numerical evidence and 
general arguments that in a ID system of lattice bosons 
with strong off-diagonal disorder the criticality of the SF- 
MG transition is the same as for the SF-MI transition in a 
pure system. We did not find the evidence in favor of the 



strong-randomness fixed point 'p| . We argued, that an al- 
ternative scenario, if in principle possible, is inconsistent 
with continuously vanishing superfiuid stiffness or/and 
compressibility at the critical point, and must predict a 
jump of the Luttinger-liquid parameter with the ampli- 
tude larger or equal to 2. 
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